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Popular Summary 


To improve our understanding of global energy and water cycle variability, and to improve 
model simulations of climate variations, it is vital to have accurate latent heat fluxes (LHF) 
over global oceans. Monthly LHF, 10-m wind speed (U IOm ), 10-m specific humidity (Q 10m ), 
and sea-air humidity difference (Q s -Q 10m ) of GSSTF2 (version 2 Goddard Satellite-based 
Surface Turbulent Fluxes) over global oceans during 1992-93 are compared with those of 
HOAPS (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data), NCEP 
(NCEP/NCAR reanalysis), and da Silva (da Silva et al.). The mean differences, standard 
deviations of differences, and temporal correlation of these monthly variables over global 
oceans during 1992-93 between GSSTF2 and each of the three datasets are analyzed. The 
large-scale patterns of the 2yr-mean fields for these variables are similar among these four 
datasets, but significant quantitative differences are found. 

The temporal correlation is higher in the northern extratropics than in the south for all 
variables, with the contrast being especially large for da Silva as a result of more missing ship 
data in the south. The da Silva has extremely low temporal correlation and large differences 
with GSSTF2 for all variables in the southern extratropics, indicating that da Silva hardly 
produces a realistic variability in these variables. The NCEP has extremely low temporal 
correlation (0.27) and large spatial variations of differences with GSSTF2 for Q s -Q 10m in the 
tropics, which causes the low correlation for LHF. Over the tropics, the HOAPS LHF is 
significantly smaller than GSSTF2 by -31% (37 W m' 2 ), whereas the other two datasets are 
comparable to GSSTF2. This is because the HOAPS has systematically smaller LHF than 
GSSTF2 in space, while the other two datasets have very large spatial variations of large 
positive and negative LHF differences with GSSTF2 to cancel and to produce smaller 
regional-mean differences. Our analyses suggest that the GSSTF2 latent heat flux, surface air 
humidity, and winds are likely to be more realistic than the other three flux datasets examined, 
although those of GSSTF2 are still subject to regional biases. 
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ABSTRACT 

Monthly latent heat fluxes (LHF), 10-m wind speed (U 10m ), 10-m specific humidity (Q ]0m ), and 
sea-air humidity difference (Q s -Q 10m ) of GSSTF2 (version 2 Goddard Satellite-based Surface 
Turbulent Fluxes) over global oceans during 1992-93 are compared with those of HOAPS 
(Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data), NCEP (NCEP/NCAR 
reanalysis), and da Silva (da Silva et al.). The mean differences, standard deviations of differences, 
and temporal correlation of these monthly variables over global oceans during 1992-93 between 
GSSTF2 and each of the three datasets are analyzed. The large-scale patterns of the 2yr-mean 
fields for these variables are similar among these four datasets, but significant quantitative 
differences are found. 

The temporal correlation is higher in the northern extratropics than in the south for all variables, 
with the contrast being especially large for da Silva as a result of more missing ship data in the 
south. The da Silva has extremely low temporal correlation and large differences with GSSTF2 for 
all variables in the southern extratropics, indicating that da Silva hardly produces a realistic 
variability in these variables. The NCEP has extremely low temporal correlation (0.27) and large 
spatial variations of differences with GSSTF2 for Q s -Q IOm in the tropics, which causes the low 
correlation for LHF. Over the tropics, the HOAPS LHF is significantly smaller than GSSTF2 by 
-31% (37 W m' 2 ), whereas the other two datasets are comparable to GSSTF2. This is because the 
HOAPS has systematically smaller LHF than GSSTF2 in space, while the other two datasets have 
very large spatial variations of large positive and negative LHF differences with GSSTF2 to cancel 
and to produce smaller regional-mean differences. Our analyses suggest that the GSSTF2 latent 
heat flux, surface air humidity, and winds are likely to be more realistic than the other three flux 
datasets examined, although those of GSSTF2 are still subject to regional biases. 
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1. Introduction 

Information on the freshwater (precipitation-evaporation), heat, and momentum fluxes at the sea 
surface is essential in understanding the interactions between the atmosphere and oceans, as well as 
global energy and water cycle variability, and in improving model simulations of climate variations. 
To predict or understand the sea surface temperature (SST) change that crucially affects the global 
climate, it is paramount that we know the variability of net surface heat flux. The spatial variations 
of seasonal and interannual variability of net surface heat flux are dominated by the variability of 
latent heat flux (LHF). This is because the spatial variation of LHF is significantly larger than that 
of shortwave radiative flux, although the magnitude of solar heating is larger than that of evaporative 
cooling (e.g., da Silva et al. 1994; Josey et al. 1999; Chou et al. 2003b). Thus it is important to 
have a reliable long-term global dataset of LHF. 

Surface measurements of LHF and other fluxes are scarce in both space and time. The 
Comprehensive Ocean-Atmosphere Data Set (COADS) has collected the most complete surface 
marine observations since 1854, mainly from merchant ships (Woodruff et al. 1993). However, the 
COADS-based fluxes have serious spatial and temporal sampling problems (measurements over 
the ship lanes with fair weather bias), together with measurement uncertainty (e.g., da Silva et al. 
1994; Josey et al. 1999; Wang and McPhaden 2001). Therefore, it is desirable that long-term 
global datasets of these fluxes be derived either from satellite observations or general circulation 
models (GCM). The World Climate Research Program (WCRP)/Global Energy and Water 
Experiment (GEWEX) Radiation Panel has established a SEA surface turbulent FLUX project, 
called SEAFLUX (Curry et al. 2003). The goal is to produce accurate high-resolution satellite- 
based datasets of surface turbulent fluxes (momentum, latent heat, and sensible heat) over global 
oceans. 

Currently, there are several datasets of global ocean LHF available, which are based on the 
surface air humidity and winds derived from the Special Sensor Microwave/Imager (SSM/I) on 
board a series of the Defense Meteorological Satellite Program (DMSP) spacecraft. The Hamburg 
Ocean Atmosphere Parameters and Fluxes from Satellite Data (HO APS; Gra(3l et al. 2000) has 
daily and monthly LHF over global oceans with 0.5° spatial resolution for the period July 
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1987-December 1998, based on the method of Schulz et al. (1997). The Goddard Satellite-based 
Surface Turbulent Fluxes (GSSTF) has two versions of global ocean LHF derived from the SSM/I 
radiances. The version 1 (GSSTF1) has daily and monthly fields for July 1987-December 1994 
with a spatial resolution of 2.0° x 2.5° lat-lon (Chou et al. 1997, 2000). The version 2 (GSSTF2) 
has daily and monthly fields for July 1987-December 2000 with 1° resolution (Chou et al. 2003). 
The Japanese Ocean Flux dataset with Use of Remote sensing Observations (J-OFURO) has 
monthly turbulent heat fluxes over global oceans with 1° resolution for 1991-95 (Kubota et al., 
2002). The J-OFURO has been further extended to the period 1992-2000 with 3 days temporal 
resolution. These flux datasets and others have been distributed to the SEAFLUX web site for 
intercomparison studies (Curry et al. 2003). 

The LHF is derived using various bulk flux algorithms from surface winds, surface air humidity 
and temperature, and SST, all of which may have a large uncertainty in the reanalyses, satellite 
retrievals, and COADS. There is no "ground truth" for the global LHF fields; thus it is important 
to conduct intercomparison studies to assess the sources of errors for various global LHF products. 
The studies can identify the strengths and weaknesses of various LHF products, and provide 
important information for improving atmospheric GCMs and satellite retrievals. As the 
NCEP/NCAR reanalysis (Kalnay et al. 1996, referred to as NCEP) and da Silva et al. (1994, 
referred to as da Silva) have been widely used for various climate studies, it is important to include 
them for the comparison. Kubota et al. (2003) compared the LHF of GSSTF1, HO APS, J- 
OFURO, NCEP, and the ECMWF (the European Centre for Medium-Range Weather Forecasts) 
analysis for 1992-94, as well as the LHF of da Silva for 1992-93, over the global oceans. They 
found that the large-scale patterns of LHF are generally similar but with significant quantitative 
differences among various products. 

Chou et al. (2003) compared the zonal averages of the GSSTF2 LHF and input parameters over 
global oceans with those of GSSTF1, HOAPS, NCEP, and da Silva for the 2yr-mean of 1992-93. 
They found that these datasets had significant differences in these parameters with GSSTF2, except 
GSSTF 1. For the GSSTF1, the differences with GSSTF2 were generally negligibly small for the 
input parameters, but not for LHF. The difference in LFEF was mainly caused by a different von 
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Karman constant for humidity between GSSTF1 (0.45) and GSSTF2 (0.40). Averaging the LHF 
over the oceans within 60°S-60°N, the GSSTF1 was 12.7 W m' 2 larger (-12% of the global 
averaged LHF of GSSTF2) than GSSTF2. The purpose of this paper is to present a detailed 
comparison study of the GSSTF2 LHF and input parameters with those of HO APS, NCEP and da 
Silva over global oceans during 1992-93. Section 2 describes the data sources used in this study. 
Section 3 briefly discusses the method for deriving GSSTF2 LHF and retrieval accuracy. Section 4 
presents the mean differences, standard deviations of differences, and temporal correlation of 
monthly input parameters for LHF (10-m wind speeds, 10-m specific humidity, and sea-air 
humidity differences) over global oceans during 1992-93 between GSSTF2 and each of the three 
datasets. The same analyses for LHF are discussed in section 5. Concluding remarks are given in 
section 6. 

2. Data sources 

The basic data used in this study are 1° x 1 0 lat-lon monthly-mean latent heat fluxes, 10-m wind 
speed (U 10m ), 10-m specific humidity (Q 10m ), and sea-air humidity difference (Q s -Qi 0m ) over global 
oceans for the period 1992-93 taken from GSSTF2, HOAPS, NCEP, and da Silva. For proper 
comparison, the 2-m specific humidity of NCEP, and 20-m wind speed and 20-m specific humidity 
of da Silva are adjusted to thelO-m height using the GSSTF2 bulk flux model. In addition, only the 
space and time matched monthly mean valid data for the common period of 1992-93 are used for 
the comparison. 

To validate GSSTF2 daily LHF and input parameters, hourly measurements of surface 
meteorology and LHF of nine field experiments conducted by the NOAA/Environmental 
Technology Laboratory (ETL) research ships over the tropical and midlatitude oceans during 
1991-99 are used (Fairall et al. 1997, 2003; Brunke et al. 2003). Table 1 shows the periods and 
locations of these nine experiments: the Atlantic Stratocumulus Transition Experiment (ASTEX), 
the Coupled Ocean-Atmosphere Response Experiment (COARE), the Fronts and Atlantic Storm 
Track Experiment (FASTEX), the Joint Air-Sea Monsoon Interaction Experiment (JASMINE), the 
Kwajalein Experiment (KWAJEX), the Nauru ’99 (NAURU99), the Tropical Instability Wave 
Experiment (TTWE), the Pan-American Climate Study in the eastern Pacific during 1999 
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(PACSF99), and the buoy service in the North Pacific (MOORINGS). The experiments provide 
hourly (50 min) covariance latent heat flux derived using the covariance or eddy correlation method. 
To reduce the flow distortion effects, only the data with the relative wind direction within 30° of the 
bow are used (Yelland et al. 1998). 

3. Derivation and validation of GSSTF2 LHF 

The GSSTF2 flux model is a bulk aerodynamics algorithm based on the Monin-Obukhov 
similarity theory including the salinity and cool-skin effects. Chou et al. (2003) have discussed this 
model somewhat in detail. The daily-mean input parameters for turbulent flux calculations include 
SSM/I U 10m of Wentz (1997) and the Q 10m , which is derived from the SSM/I precipitable water of 
the entire atmospheric column and the 500-m layer near the surface (Chou et al. 1995, 1997). The 
input parameters also include the 2-m air temperature and SST of the NCEP/NCAR reanalysis 
(Kalney et al. 1996). The saturation specific humidity at the sea surface (Q s ) for computing LHF 
depends on the skin SST. The skin SST is generally cooler than the bulk SST with a daily mean 
difference of ~0.2 °C, but its instantaneous value relies on various bulk-skin SST difference 
models, surface net heat flux, and wind-induced oceanic mixing (e.g. Fairall et al. 1996b; Webster 
et al. 1996; Wick et al. 1996; Curry et al. 2003). In addition, the saturation vapor pressure is - 2 % 
smaller for saline water than for pure water (Fairall et al. 1996a; Zeng et al. 1998). To partially 
compensate for the salinity and cool skin effects on LHF, the GSSTF2 flux model estimates Q s 
using the approximated formulation as Q s = 0.98 (0.622 e/P). Here e is the saturation vapor 
pressure for pure water at the bulk SST, and P is the sea level pressure taken from the 
NCEP/NCAR reanalysis. The GSSTF2 neutral moisture transfer coefficient at 10-m height is 
1.05-1.22 x 10 3 for U ]0m of 3-18 m s' 1 . The coefficient increases with decreasing winds for U 10m 
< 3 m s' 1 . The GSSTF2 moisture transfer coefficient (for U 10m up to -18 m s' 1 ) is in close 
agreement with those of Zeng et al. (1998), Fairall et al. (1996a), and Renfrew et al. (2002). 

Chou et al. (2003) validated the GSSTF2 bulk flux model, using hourly input parameters from 
ship measurements of ten field experiments conducted by the NOAA/ETL. They found that the 
hourly LHF computed using the GSSTF2 flux model agreed very well with observations. 
Compared to 1913 samples of hourly observed LHF of the ten experiments, the computed hourly 
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LHF had a bias of 4.5 W m' 2 , a standard deviation (SD) error of 19.6 W m' 2 , and a correlation of 
0.91, with a correlation range of 0.77-0.95. The SD error is the standard deviation of the 
differences between the computed and observed. Their results suggested that the GSSTF2 flux 
model is generally accurate for weak and moderate winds, but slightly underestimates LHF for 
strong winds. 

Chou et al. (2003) also validated the GSSTF2 daily turbulent fluxes and input parameters with 
the nine NOAA/ETL field experiments shown in Table 1. The comparison of daily LHF and 
relevant input parameters between GSSTF2 and the nine field experiments are shown in Table 2 
and Fig. 1. Using the GSSTF2 flux model, the GSSTF2 daily surface wind speed (U) and surface 
air specific humidity (Q a ) are adjusted to the measurement heights of the ships (~ 14-21 m) for 
proper validation. The GSSTF2 daily SST is compared with the bulk SST at the 5-cm depth 
measured by the ships. The comparison with five major tropical experiments with large samples 
(ASTEX, CO ARE, JASMINE, NAURU99, and KWAJEX) is also shown in Table 2 for reference. 

Compared to the nine (five tropical) experiments. Fig. 1 and Table 2 show that, the GSSTF2 
daily LHF has a bias of 0.8 (-2.6) W m' 2 , a SD error of 35.8 (29.7) W m' 2 , and a correlation of 
0.83 (0.80). The GSSTF2 daily wind speed has a bias of 0.36 (0.31) m s' 1 , a SD error of 1.38 
(1.07) m s' 1 , and a correlation of 0.92 (0.87). The GSSTF2 daily surface air specific humidity has 
a bias of 0.67 (1.01) g kg' 1 , a SD error of 1.23 (1.11) g kg' 1 , and a correlation of 0.97 (0.85). The 
GSSTF2 daily SST has a bias of 0.04 (0.02) °C, a SD error of 0.51 (0.30) °C, and a corrlelation of 
1.0 (0.99). We note that the GSSTF2 daily Q a has a positive bias of ~1 g kg' 1 in the tropical 
oceans, especially for Q a >16 g kg' 1 (Fig. lc). However, the LHF has a rather small negative bias of 
— 2.6 W m' 2 in the tropical oceans. This is mainly due to the fact that, in the tropical oceans, the 
GSSTF2 flux model produces a small positive bias of ~5 W m' 2 for LHF (Chou et al. 2003) and 
that the weak winds associated with moist air have a small positive bias of ~0.3 m s' 1 (Table 2). 
These two factors offset the underestimation of LHF due to a positive bias of humidity. 

Sources of retrieval-ship differences in daily LHF and input parameters include the spatial- 
temporal mismatch between GSSTF2 and ships, as well as the errors in the input parameters and 
fluxes for both GSSTF2 and ship observations. The collocated daily GSSTF2 variables are 
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computed from 2-3 satellite observations averaged over a 1° area that encloses the ship locations, 
while those of the ships are computed from at least two hourly measurements over a much smaller 
area. Assuming daily retrieval errors arc independent. Table 2 shows the SD errors for the 
monthly-mean LHF reduce to 6.5 and 5.4 W m' 2 , as inferred from the nine and five tropical 
experiments, respectively. The SD errors for the monthly-mean U, Q a , and SST reduce to 0.25 
(0.20) m s' 1 , 0.22 (0.20) g kg' 1 , and 0.09°C (0.05°C), respectively, as inferred from the nine (five 
tropical) experiments. 

4. Comparison of relevant input parameters 

a. Surface wind speeds 

The U 10m averaged over 1992-93 for GSSTF2, and the differences U 10m (HOAPS - GSSTF2), 
U 10ra (NCEP - GSSTF2), and U 10m (da Silva - GSSTF2), respectively, for the 2-yr period over 
global oceans are shown in Fig. 2. Figure 2a shows that the maximum annual mean wind speed is 
located in the trade-wind belts (~8-9 m s' 1 ) and extratropical storm-track regions (~8-l 1 m s' 1 ). 
The minimum wind speed is located in the weak-wind (~4-6 m s' 1 ) areas of equatorial Indian 
Ocean, South Pacific convergence zone (SPCZ), and intertropical convergence zone (TTCZ), and the 
subtropical highs (~7 m s' 1 ). Figure 2b shows that the HOAPS U 10m is significantly smaller than 
GSSTF2 in the tropical oceans with a difference up to 2 m s' 1 centered around the locations of the 
maximum wind in the trade wind belts of both hemispheres. Over the extratropical oceans, on the 
contrary, the HOAPS U 10m is larger than GSSTF2 with a maximum difference of 2 and 1.5 m s' 1 , 
centered near 50°N and 50°S, respectively. 

Figure 2c shows that the NCEP U 10m is also significantly smaller than GSSTF2 in the tropics 
(equatorward of -40°, except the eastern boundary current regions in the Pacific and Atlantic), with 
the maximum difference of 2 m s' 1 near the equator. Over the extratropics, the NCEP U 10m is larger 
than GSSTF2 with the maximum difference of 0.5 and 1 m s' 1 around 50°N and 50°S, respectively. 
Figure 2d shows that Ui 0m (da Silva - GSSTF2) has small-scale features and is quite different with 
the large-scale patterns of U 10m (HOAPS - GSSTF2) and U 10m (NCEP - GSSTF2) shown in Figs. 
2b and 2c. Over the oceans north of ~10°N, U 10m (da Silva - GSSTF2) is generally positive in the 
high wind areas of extratropical storm track regions and trade wind belts (by ~2-3 m s' 1 ), but is 
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generally negative in the weak wind areas of subtropical highs (by ~l-2 m s' 1 ). Over the oceans 
south of 10°N, it characterizes the small-scale features with neighboring significant large positive 
and negative difference centers, especially with the magnitude of U 10m difference reaching ~4-5 m s' 
1 over the data void region of the eastern South Pacific. 

Figure 3 shows the standard deviations of differences (SDD) and temporal cross correlation of 
monthly U 10m between GSSTF2 and each of HOAPS, NCEP, and da Silva, respectively, during 
1992-93. The SDD fields are shown in the left column, while those of correlation are shown in the 
right column. Small SDD indicates a small temporal variability for the monthly differences 
between datasets, or small dispersions from the mean difference fields shown in Fig. 2. Figures 3a 
and 3b show that, over most of the tropical oceans, the SDD in U 10m generally has a small value of 
-0.5-0.7 m s' 1 with high correlation of >0.9, which implies that monthly U 10m (HOAPS - GSSTF2) 
has a small temporal variability there. This indicates that the U 10m of HOAPS and GSSTF2 have a 
large systematical difference (Fig. 2b) but with similar temporal variability in the tropical oceans. In 
the extratropical oceans, the SDD in U IOm is generally -0.5-1. 0 m s' 1 with high correlation of 
-0.8-0.9 (-0.8) in the north (south), which implies that monthly U 10m (HOAPS - GSSTF2) has a 
slightly increased temporal variability there. 

Figures 3c and 3d show that, poleward of -10°, the SDD in U 10m generally has a small value of 
-0.5-0.7 m s' 1 with high correlation of >0.9 (0.8-0.9) in the north (south), which implies that 
monthly U 10m (NCEP - GSSTF2) has a small temporal variability there during 1992-93. Within 
-10° of equator, the SDD in U, 0m is generally -0.7-1 .0 m s' 1 with the correlation reduced to 
-0.3-0.6, which implies that monthly U 10m (NCEP - GSSTF2) has a large temporal variability 
there. Figures 3e and 3f show that the SDD in U 10m generally has a large value of -1-1.5 m s' 1 
over most of the global oceans with poor correlation of <0.6 in the oceans south of 30°N (-0. 8-0.9 
in the ocean north of 30°N)- This indicates that the monthly U 10m of da Silva and GSSTF2 do not 
correlate well and that monthly U 10m (da Silva - GSSTF2) has a large temporal variability over most 
of the global oceans, especially over the large data-void areas of the Southern Hemisphere. 

The Wentz (1997) SSM/I U, 0m from 1987 to 1997 have been extensively evaluated with those 
of the Tropical Atmosphere-Ocean (TAO) and National Data Buoy Center (NDBC) buoys by 
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Mears et al. (2001), and those of the ECMWF analysis and NCEP by Meissner et al. (2001). 
Meissner et al. (2001) pointed out that both global analyses did not assimilate the Wentz wind 
products and that the SSM/I wind speeds assimilated in NCEP were derived using a neural network 
algorithm and were different from those of Wentz (1997). Mears et al. (2001) found that the mean 
difference between SSM/I and buoy winds was typically < 0.4 m s' 1 and the SD error was <1.4 m s' 
*. Meissner et al. (2001) found that the collocated SSM/I and NCEP U 10m had a SD difference of 
2.4 m s' 1 , which reduced to 1.2 m s' 1 for the monthly averages, and that the NCEP U 10m was 
underestimated in the tropical Pacific and tropical Atlantic. Wang and McPhaden (2001) found the 
NCEP surface winds were weaker than those of TAO buoys by -1-1.5 m s' 1 in the tropical Pacific. 
Smith et al. (2001) found the NCEP U 10m was underestimated by 0.4-1. 0 m s' 1 for five geographic 
regions (the North Atlantic, South Atlantic, Southern Ocean, Western Pacific, and Arabian Sea), as 
compared to those measured by research ships during the 1990-95 World Ocean Circulation 
Experiment. Renfrew et al. (2002) found that the NCEP U 10m was 0.4 m s' 1 higher than those 
measured by the research ship over the Labrador Sea. These results are consistent with this study 
(Table 2 and Fig. 2c). 

Da Silva et al. (1994) assumed an anemometer height of 20 m to the entire wind dataset 
measured by ship anemometers of CO ADS to derive a Beaufort equivalent scale for determining 
visual wind speeds, which depend on sea states. However, Kent and Taylor (1997) found that the 
true anemometer heights had large standard deviations with the means generally much higher than 
20 m and increasing with time. For example, they found that the mean anemometer height was 35.2 
m (24.2 m) with a standard deviation of 8.4 m (10.9 m) in middle latitudes of the North Pacific 
(North Atlantic) during 1990. An underestimation of anemometer height can cause unrealistic 
higher ship anemometer-measured (and visual) wind speeds, because the higher wind speeds 
measured at the higher anemometer heights are assigned to the assumed lower levels. For the same 
error of anemometer height, the stronger the wind is, the larger the wind speed error is. This can 
cause a larger discrepancy of wind speeds in the high-wind regions. This is likely to be the major 
reason that a large discrepancy of U 10m in the high-wind region of the Northern Hemisphere is 
found between da Silva and GSSTF2. 


Table 2 and Fig. lb suggested that the GSSTF2 monthly wind speed has a bias of 0.36 (0.31) 
m s' 1 , and a SD error of 0.25 (0.20) m s' 1 , respectively, as inferred from the nine (five tropical) 
experiments. This result and the above mentioned previous studies suggest that the GSSTF2 U 10m 
is more accurate and that the differences and SDD of U IOm shown in Figs. 2 and 3 are most likely 
mainly caused by the errors in HO APS, NCEP, and da Silva, 
b. Surface air specific humidity 

Figure 4 shows the Q 10m averaged over 1992-93 for GSSTF2, and the differences 
Q 10m (HOAPS - GSSTF2), Q 10m (NCEP - GSSTF2), and Q 10m (da Silva - GSSTF2), respectively, 
over global oceans. Figure 4a shows that the Q 10m has a maximum of >18 g kg' 1 in the equatorial 
warm pool and decreases poleward, with the pattern similar to that of SST. The large-scale patterns 
of Q 10m for HO APS, NCEP, and da Silva averaged over 1992-93 (not shown) are similar to that of 
GSSTF2. However, Figs. 4b-4d show significant differences with GSSTF2. Figure 4b shows that 
the HO APS Q 10nl is higher than GSSTF2 over the global oceans with the maximum difference of 
1. 0-1.5 g kg' 1 in the tropics with high SST. The minimum Q 10m difference of <~0.5 g kg' 1 is 
located in the areas poleward of -50°, as well as the dry tongue regions of the eastern South Pacific 
and Atlantic. 

Figure 4c shows that, over the equatorial Indian Ocean, SPCZ, and ITCZ, the NCEP Q 10m is 
smaller than GSSTF2 by -0.5-1 .0 g kg' 1 , while it is significantly larger than GSSTF2 by -0.5-2.5 
g kg' 1 for the rest of the global oceans. The maximum difference in Q 10m is -2.5 g kg' 1 in the trade 
wind and dry tongue area of the eastern South Pacific, -2 g kg' 1 in the other trade wind areas and 
-1 g kg' 1 in the extratropical oceans. We note that the dry tongue of the HOAPS Q 10m is similar to 
that of GSSTF2 over the eastern South Pacific, but the dry tongue of NCEP is significantly weaker 
than that of GSSTF2 (not shown). This result can be inferred from Figs. 4b and 4c. Figure 4d 
shows that Q 10m (da Silva - GSSTF2) is characterized by small-scale features with large magnitudes 
of differences over data sparse regions, especially the equatorial Pacific and oceans south of ~30°S. 
The large-scale pattern of Q IOm (da Silva - GSSTF2) (Fig. 4d) is similar to that of Q 10m (NCEP - 
GSSTF2) (Fig. 4c) north of ~25°S but is quite different south of 25°S. This is because the same 
observed surface air humidity data of CO ADS were used in both NCEP and da Silva. 
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We note that the Q 10m of da Silva undulates in the zonal direction (which is very different from 
the large-scale and smoothing patterns of other three data sets) and does not follow the SST pattern 
south of ~30°S (not shown). Thus there are significant smaller scale differences in the zonal belt of 
30°S-50°S (Fig. 4d). In addition, the Q 10m of da Silva is generally lower than GSSTF2 in the 
equatorial areas (where the Q 10m is the highest) but is generally higher than GSSTF2 poleward of 
these regions. Compared to GSSTF2, the da Silva Q 10m thus has smaller poleward gradient, which 
is an indication of a smoothing effect from the successive corrections of da Silva. The positive 
differences are especially large south of 30°S. Since data are scarce in the oceans south of 30°S, 
the high Q IOm (and SST) observed further equatorward is extrapolated poleward (Kent et al. 2000). 
This can cause da Silva to have unrealistic high values of Q 10m (and SST) south of 30°S and to have 
large differences with GSSTF2 there (Figs. 4d and 5d). 

Since there is no “ground truth” of surface air humidity over global oceans, we analyze the 
humidity discrepancies based on Table 2, Fig. lc, and previous studies. Table 2 shows that the 
GSSTF2 Q a has a positive bias of ~1 (~0.7) g kg' 1 in the tropical oceans (for the nine experiments 
over the tropical and midlatitude oceans). Figure lc shows that the GSSTF2 Q a has a positive bias 
for the moist region with Q a of -16-20 g kg' 1 but has a small negative bias for the dry region with 
Q a of -3-6 g kg' 1 . Thus Table 2, and Figs, lc, 4a, and Fig. 4b suggest that the HOAPS Q 10m is 
significantly overestimated in the tropical oceans, as it is larger than the moist biased Q a of 
GSSTF2. In addition, Wang and McPhaden (2001), Smith et al. (2001), and Renfrew et al. (2002) 
found that the NCEP surface air humidity had positive biases when compared with those measured 
by TAO buoys and research ships. Their results and this study imply that the positive Q 10ro (NCEP 
- GSSTF2) shown in Fig. 4c is mostly a result of the moist bias of NCEP, while the negative 
Q ]0m (NCEP - GSSTF2) is likely mostly due to the moist bias of GSSTF2. 

Previous studies (e.g., da Silva et al. 1994; Chou et al. 1997; Josey et al. 1999) found that ship 

observations overestimated dew point temperatures (by -0.5 C), which resulted in moist bias of the 
surface air humidity for CO ADS, and thus for da Silva and NCEP. Their results and this study 
suggest that the moist bias of da Silva is likely responsible for the positive Q 10m (da Silva - 


GSSTF2) shown in Fig. 4d. These analyses suggest that the GSSTF2 Q 10m is likely to have better 
quality than other three datasets analyzed, although it is still subject to regional biases. More 
extensive measurements of high quality surface meteorology over the global oceans are essential for 
further validations, 
c. Sea-air humidity differences 

Figure 5 shows the Q s -Qi 0 m averaged over 1992-93 for GSSTF2, and the differences 
[(Qs-Q.oJ(HOAPS - GSSTF2)], [(Q s -Q 10m )(NCEP - GSSTF2)], and [(Q s -Q 10 J(da Silva - 
GSSTF2)], respectively, over global oceans. It is noted that there is a 2 % reduction in Q s due to 
salinity effect in GSSTF2, HOAPS, and da Silva, but not in NCEP. In addition, there is a sort of 
cool skin effect on the GSSTF2 Q s due to the use of approximate humidity equation mentioned in 
section 3. Figure 5a shows that the Q s -Q 10m has the maximum values of ~6 g kg' 1 over the trade 
wind belts of both hemispheres and decreases to <~2 g kg' 1 poleward of ~50° resulting from the 
reduction of SST. Over the equatorial region (10°S-10°N), the Q s -Qio m has a large value of ~5 g 
kg' 1 in the high SST regions of the Indian Ocean and western Pacific, and decreases eastward in the 
eastern Pacific due to upwelling-induced decrease in the SST. Figure 5b shows that the HOAPS 
Q s -Qiom is smaller than GSSTF2 by -0.5-1 and 0-0.5 g kg' 1 , respectively, in the areas 
equatorward and poleward of -20°. The difference in Q s -Q 10m is smaller over the dry or cold 
tongue areas of the eastern South Pacific and South Atlantic, with the minimum difference ~30°S. 

Figure 5c shows that, over the regions of equatorial Indian Ocean, SPCZ, and ITCZ, the NCEP 
Q s — Q 10m is significantly larger than GSSTF2 by -0.5-1 .5 g kg' 1 . For the rest of the global oceans, 
it is significantly smaller than GSSTF2 by -0.5-2.0 g kg' 1 , with the maximum difference reaching 
-2 g kg' 1 over the trade wind area in the eastern South Pacific, up to -1.5 g kg' 1 in the other trade 
wind areas, and -1 g kg' 1 in the extratropical oceans. Figure 5d shows that the [(Q s -Q 10m )(da Silva 
- GSSTF2)] has small-scale structures with large magnitudes of differences over data sparse 
regions, especially the equatorial Pacific and the oceans south of 30°S. The large-scale pattern of 
[(Qs-Q.omXda Silva - GSSTF2)] (Fig. 5d) is similar to that of [(Q s -Q 10m )(NCEP - GSSTF2)] 
(Fig. 5c) north of ~25°S but is quite different south of 25°S. It can be seen from Figs. 4 and 5 that, 
for these datasets, the differences in Q 10m mainly control the differences in Q s -Q 10m and the 
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differences in Q s is relatively small (see Table 3), which is consistent with the results of Chou et al. 
(2003). 

Figure 6 shows the SDD and temporal cross correlation of monthly Q s -Q, 0m between GSSTF2 
and each of HOAPS, NCEP, and da Silva, respectively, during 1992-93. Figures 6a and 6b show 
that, over most of the oceans, the SDD in Q s -Q 10m is -0.4-0.7 g kg' 1 (<0.4 g kg' 1 or smaller near the 
dry tongues and polar regions) with a moderate high correlation of -0.7-0.9 (slightly smaller than 
the correlation of U 10m ). This implies that the monthly [(Q s -Q 10m )(HOAPS - GSSTF2)] has a 
moderate temporal variability over most of the oceans. Over the tropical oceans where the HOAPS 
Qs — Qiom is significantly smaller than GSSTF2 (Fig. 5b), there is high correlation associated with 
relatively small SDD in Q s -Qi 0m - This indicates that the monthly Q s -Q 10m of HOAPS is 
systematically smaller than that of GSSTF2, but with the similar temporal variability as that of 
GSSTF2. 

Figures 6c and 6d show that, over the tropical oceans, the SDD in Q s -Q 10m has a large value of 
-1-1.5 g kg' 1 with low correlation of <0.3 (reducing to about -0.3 in the equatorial Pacific and 
Indian oceans). This implies that the monthly [(Q s -Q 10m )(NCEP - GSSTF2)] has a very large 
temporal variability over the tropical oceans during 1992-93. In the extratropical oceans, the SDD 
in Q s -Q 10m decreases to -0.4-0.7 g kg' 1 and the correlation increases to -0.6-0.9. This implies that 
the monthly [(Q s -Q ]0m )(NCEP - GSSTF2)] has a reduced temporal variability over the extratropical 
oceans, as compared to over the tropical oceans. Figures 6e and 6f show that the SDD in Q s -Q 10m 
has a large value of -0.7-1. 5 g kg' 1 equatorward of ~40°N (-0.5 g kg' 1 poleward of ~40°N) with 
very low correlation of <0.3 over most of the oceans south of ~30°N. This implies that the monthly 
[(Q s — Qjom)(da Silva - GSSTF2)] has a very large temporal variability over most of the global 
oceans. Note that the low correlation associated with relatively large SDD in Q s -Q 10m for da Silva, 
as well as NCEP in the equatorial area, is generally located in the data sparse regions (e.g., da Silva 
et al. 1994; Chou et al. 1997; Josey et al. 1999). This may be due to the errors of Q 10m arising from 
the interpolation scheme of da Silva and the shortcoming of the boundary layer and cumulus 
parameterization used in NCEP, as the correlation of Q s -Q 10m (or Q IOm ) between GSSTF2 and 
HOAPS is generally high over these areas. 
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5. Comparison of LHF 

The LHF averaged over 1992-93 for GSSTF2, and the differences LHF(HOAPS - GSSTF2), 
LHF(NCEP - GSSTF2), and LHF(da Silva - GSSTF2), respectively, for the 2-yr period over 
global oceans are shown in Fig. 7. Figure 7a shows that the maximum LHF is located in the trade 
wind belts of both hemispheres (-160-180 W m' 2 ) and in the western boundary current regions of 
the Kuroshio and Gulf Stream (-160 W m' 2 ). This results from high winds (-8-9 m s' 1 ) coupling 
with large Q s -Q 10m (-5-6 g kg' 1 ) in these areas (Figs. 2a and 5a). The minimum LHF (<~60 W m‘ 
2 ) is found in the eastern equatorial Pacific and Atlantic, due to upwelling-induced cold SSTs 
associated with weak winds, and in the high latitudes due to poleward decrease of SST. We note 
that the general patterns of U 10m , Q s -Q 10m , and LHF of GSSTF2 for 1992-93 (Figs. 2a, 5a, and 7a) 
are similar to those of Chou et al. (2003) for 1988-2000. For the 1992-93 mean condition, the 
large-scale pattern of the GSSTF2 LHF shown in Fig. 7a is similar to those of HOAPS, NCEP, 
and da Silva (not shown). However, there are significant differences with GSSTF2 geographically. 

The moisture transfer coefficients of HOAPS and da Silva are derived based on Smith (1988) 
and Large and Pond (1982), respectively. Both coefficients are very close to that of GSSTF2. 
Thus, the discrepancy in the LHF among GSSTF2, HOAPS, and da Silva is primarily caused by 
the differences in the input parameters. However, this situation is not true for NCEP. The moisture 
transfer coefficient of NCEP is significantly larger than that of GSSTF2 (and the other two flux 
datasets; Zeng et al. 1998), which appears to offset the effects of weaker winds and smaller sea-air 
humidity difference on the LHF (Figs. 2c and 5c). 

Figure 7b shows that the HOAPS LHF is significantly smaller than that of GSSTF2 over the 
global oceans, especially in the tropics. In the tropical oceans, the difference is -20-60 W m' 2 , with 
the maximum differences of -60 and 40 W m' 2 centered ~10°N and 10°S, respectively. The 
difference decreases poleward and is negligible south of ~30°S. In the tropical oceans, the negative 
LHF(HOAPS - GSSTF2) is a result of smaller U IOm and Q s -Q 10m of HOAPS (Figs. 2b and 5b). 
For the negative LHF(HOAPS - GSSTF2) in the oceans north of 30°N, the effect of the smaller 
HOAPS Q s -Q 10m (coupling with high winds) slightly overcompensates the effect of higher U 10m 
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(coupling with low Q s -Q 10m ) on the LHF. For the negligible LHF difference in the oceans south of 
~30°S, the effects of smaller Q s -Q 10m and larger U 10ra on LHF generally offset each other. 

Figure 7c shows that the NCEP LHF is larger than GSSTF2 in the equatorial Indian Ocean, 
SPCZ, ITCZ, the Kuroshio and Gulf Stream, and subtropics, with the maximum difference up to 40 
W m' 2 . For the rest of the global oceans, the NCEP LHF is smaller than GSSTF2, with the 
maximum difference of ~60 W m‘ 2 located in the dry tongue and trade wind region of the eastern 
South Pacific. For the trade wind regions (especially over the South Pacific dry tongue), the 
negative LHF(NCEP - GSSTF2) is a result of smaller U 10m and Q s -Q 10m of NCEP (Figs. 2c and 
5c). The positive LHF(NCEP - GSSTF2) in the equatorial Indian Ocean, SPCZ, and ITCZ is due 
to the effects of larger NCEP Q s -Q 10m and moisture transfer coefficient overcompensating the 
effect of smaller U 10m on the LHF. The negative LHF(NCEP - GSSTF2) in the extratropical 
oceans is due to the effect of smaller NCEP Q s -Qi 0m overcompensating the effects of larger U 10m 
and moisture transfer coefficient on the LHF. 

Figure 7d shows the LHF(da Silva - GSSTF2) has small-scale structures, with its large-scale 
pattern somewhat similar to that of LHF(NCEP - GSSTF2) (Fig. 7c). For example, the da Silva 
LHF is larger than GSSTF2 in the equatorial Indian Ocean, SPCZ, ITCZ, and the Kuroshio area, 
but with a larger difference of 20-60 W m' 2 . In addition, the da Silva LHF is smaller than GSSTF2 
in the trade wind regions and the oceans south of 30°S, but with significantly larger localized 
differences of 20-80 W m' 2 . Note that Fig. 7d shows that there are some neighboring extreme 
large positive and negative difference centers of LHF located in the data sparse regions for da Silva. 
This feature is quite different than the more organized large-scale patterns of LHF difference for 
HO APS and NCEP (Figs. 7b and 7c). 

Figure 8 shows the SDD and temporal cross correlation of monthly LHF between GSSTF2 
and each of HO APS, NCEP, and da Silva, respectively, during 1992-93 over global oceans. Figure 
8a and 8b show that the SDD in LHF is generally -10-20 W m' 2 over the oceans, with a high 
correlation of 0.8-0.9 (0.6-0.8) in the oceans north of ~20°S (south of 20°S). This implies that the 
monthly LHF(HOAPS - GSSTF2) has a relatively smaller temporal variability in the oceans north 
of 20°S than the ocean south of 20°S. 


Figure 8c and 8d show that, over the tropical oceans, the SDD in LHF has a large value of 
-20-30 W m" 2 with a low correlation of 0.6-0.9 (0.3-0.6 in the ITCZ), which implies that the 
monthly LHF (NCEP - GSSTF2) has a large temporal variability there, especially over the ITCZ. 
In the extratropical oceans, the SDD in LHF reduces to -10-20 W m' 2 with an increased 
correlation of 0.8-0.9 (larger for the northern than the southern extratropical oceans). This implies 
that the monthly LHF (NCEP - GSSTF2) has a smaller temporal variability in the extratropical 
oceans than in the tropical oceans. Figure 8e and 8f show that the SDD in LHF generally has a 
large value of -20-40 W m' 2 over most of the global oceans, with the correlation of <0.6 over most 
of the areas south of ~30°N. This implies that the monthly LHF(da Silva - GSSTF2) has a large 
temporal variability over most of the oceans, especially the equatorial and southern extratropical 
oceans 

Table 3 shows the regional averages of LHF and input parameters for the global oceans 
(60°S-60°N), northern extratropical oceans (20°N-60°N), tropical oceans (20°S-20°N), and 
southern extratropical oceans (20°S-60°S) during 1992-93 for GSSTF2, as well as the differences 
of HO APS, NCEP, and da Silva from GSSTF2. The global-mean LHF is the largest for GSSTF2 
(108.2 W m' 2 ) and is the smallest for HO APS (88.5 W m' 2 ), with a difference of 20 W m' 2 . Over 
the tropical oceans, the HO APS LHF is significantly smaller than GSSTF2 by -31% (37 W m' 2 ), 
whereas the other two datasets are comparable to GSSTF2. This is because the HO APS LHF is 
systematically smaller than GSSTF2 in space, while the other two datasets have veiy large spatial 
variations in LHF with large positive and negative differences, which cancel to produce smaller 
regional-mean differences (Fig. 7). The smallness of HO APS LHF in the tropical oceans is mainly 
a result of weaker U 10m (by -1.1 m s' 1 ) coupling with smaller Q s -Q 10m (by -0.7 g kg' 1 ). The 
significant larger Q 10m (by 1.1 g kg' 1 ) is the main cause for the smallness of HO APS Q s -Qi 0m - 

The global-mean LHF of NCEP is comparable to that of GSSTF2, but Ui 0m and Q s -Q 10m are 
smaller (by -0.6 m s' 1 and -0.5 g kg' 1 ), which appear to offset the larger moisture transfer 
coefficient (Zeng et al. 1998). Due to the cancellation of large local positive and negative 
differences with GSSTF2, the regional-mean LHF of da Silva appears to be comparable to that of 
GSSTF2, except for the southern extratropical oceans. Over the southern extratropics, the da Silva 



LHF is the smallest among the four datasets (Fig. 7 and Table 3), -23% (22 W m' 2 ) smaller than 
GSSTF2. This discrepancy is mainly due to the larger Q 10m (by -1.2 g kg' 1 ) and is most likely due 
to the errors arising from the interpolation of missing data in the large data-void regions (from the 
relatively high Q 10m in the low to mid-latitudes) as discussed before (see Fig. 4d). 

Table 4 shows regional-mean standard deviation of differences/temporal cross correlation of 
monthly LHF and input parameters for the global oceans, northern extratropical oceans, tropical 
oceans, and southern extratropical oceans between GSSTF2 and each of HOAPS, NCEP, and da 
Silva during 1992-93. Note that, for the same region, a low SDD value is associated with a high 
correlation, indicating similar behavior of two time series. Averaged over the global oceans, 
HOAPS has the highest correlation in LHF with GSSTF2 (0.81), due to the highest correlation in 
U 10m and especially Q s -Q 10m - On the other hand, da Silva has the lowest correlation in LHF with 
GSSTF2 (0.51), due to the lowest correlation in U 10m and especially Q s -Q, 0m - The correlation is 
higher in the northern extratropical oceans than in the south for all variables, with the contrasts 
especially large for da Silva as a result of more missing data in the south. Over the tropical oceans, 
NCEP has a very low correlation in LHF with GSSTF2 (0.69). This is primarily a result of very 
low correlation in Q s -Q 10m (0.27), which is close to that of da Silva (0.31) but is significandy 
smaller than that of HOAPS (0.81). This may be due to missing data and the shortcoming of the 
cumulus and boundary parameterization used in NCEP. In short, our analyses suggest that the 
GSSTF2 LHF, surface air humidity, and winds are likely to be more realistic than the other three 
datasets analyzed, although those of GSSTF2 are still subject to regional biases. More high-quality 
observations over global oceans are needed to do a more detailed regional validation and to confirm 
our conclusion. 

6. Concluding remarks 

The monthly LHF and input parameters (U 10m , Q 10m , and Q s -Qi 0m ) of GSSTF2 over global 
oceans during 1992-93 are compared with those of HOAPS, NCEP, and da Silva. The large-scale 
patterns of the 2yr-mean fields for these variables are similar among these four datasets, but 
significant quantitative differences are found. Over the tropics, the HOAPS LHF is significantly 
smaller than GSSTF2 by -31% (37 W m' 2 ), whereas the other two datasets are comparable to 
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GSSTF2. This is because that HOAPS has systematically smaller LHF than GSSTF2 in space, 
while the other two datasets have very large spatial variations of large positive and negative LHF 
differences with GSSTF2 which cancel to produce smaller regional-mean differences. The 
smallness of the HOAPS LHF is a result of smaller Q s -Q 10m and U 10m , with the former mainly due 
to a larger Q 10m (by- 1.1 g kg 1 ). Over the northern (southern) extratropical oceans, the HOAPS 
LHF is smaller than GSSTF2 by 14 W m' 2 (2 W m' 2 ). Averaged over the global oceans 
(60°S-60°N), the HOAPS LHF is smaller than GSSTF2 by 20 W m' 2 . 

In the equatorial Indian Ocean, SPCZ, ITCZ, the Kuroshio and Gulf Stream, and subtropics, the 
NCEP LHF is larger than GSSTF2, with the maximum difference up to 40 W m' 2 . For the rest of 
the oceans, it is smaller than GSSTF2, with the maximum difference of -60 W m' 2 located in the 
dry tongue and trade wind region of the eastern South Pacific. The smaller NCEP LHF in the trade 
wind regions is a result of smaller Q s -Q 10 m an ^ U 10m . For the equatorial Indian Ocean, SPCZ, and 
ITCZ, the larger NCEP LHF is because the effects of larger NCEP Q s -Q 10m and moisture transfer 
coefficient overcompensate the effect of smaller U 10m on the LHF. The smaller NCEP LHF in the 
extratropical oceans is because the effect of smaller NCEP Q s -Q 10m overcompensates the effects of 
larger U 10m and moisture transfer coefficient on the LHF. Averaged over the global oceans, the 
NCEP LHF is smaller than GSSTF2 by only -3 W m' 2 due to the large cancellation of local 
positive and negative differences with GSSTF2. The NCEP has extremely low temporal correlation 
(0.27) and large spatial variations of differences with GSSTF2 for Q s -Q 10m in the tropics, which 
causes the low correlation for LHF. The very low temporal correlation of Q s -Q IOm for NCEP, 
which is generally located in the data sparse tropics, is close to that of da Silva (0.31) but is 
significantly smaller than that of HOAPS (0.81). 

The LHF difference of da Silva with GSSTF2 has small-scale structures, with its large-scale 
pattern somewhat similar to that of NCEP, but with larger magnitudes of difference. The da Silva is 
also found to have some nearby extreme large positive and negative LHF difference centers located 
in the data sparse regions. This feature is quite different than the more organized large-scale LHF 
difference patterns of HOAPS and NCEP (with GSSTF2). Due to the cancellation of these large 
local positive and negative differences, the regional-mean LEIF of da Silva appears to be comparable 
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to that of GSSTF2, except for the southern extratropical oceans with -23% (22 W m' 2 ) smaller 
value than GSSTF2. This discrepancy is mainly due to a larger Q 10m (by -1.2 g kg' 1 ) and is most 
likely due to the errors arising from the interpolation in the large data-void regions. In addition, da 
Silva has extremely low temporal correlation and large differences with GSSTF2 for all variables in 
the southern extratropics, indicating that da Silva hardly produces a realistic variability in these 
variables. The temporal correlation is higher in the northern extratropics than in the south for all 
variables, with the contrast being especially large for da Silva as a result of more missing ship data 
in the south. 

Our analyses suggest that the GSSTF2 LHF, surface air humidity, and winds are likely to be 
more realistic than the other three flux products examined, although those of GSSTF2 are still 
subject to regional biases. This conclusion is consistent with that of Chou et al. (2003). More 
high-quality observations covering wider areas of the global oceans (or more areas covering 
different climatic regimes) are vital to do a more complete regional validation, to improve satellite 
retrieval, and to further confirm our conclusion. The GSSTF2 LHF is useful for climate studies 
and is available at http://daac.gsfc.nasa.gOv/CAMPAIGN_DOCS/hydrology/hd_gsstf2.0.html. 
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Table 1. Times and locations of nine field experiments conducted by the NOAA/ETL research 
ships. 


Experiments 

Times 

Locations 

ASTEX 

92/6/6-92/6/28 

30°N, 36°W 

COARE 

92/11/11-93/2/16 

1.7°S, 156°E 

FASTEX 

96/12/22-97/1/26 

42-52°N, 5-60°W 

JASMINE 

99/5/4-99/5/31 

5°S-13°N, 88-98°E 

KWAJEX 

99/7/28-99/9/10 

9°N, 167°E 

MOORINGS 

99/9/14-99/10/21 

8°N, 167°E-49°N, 130°W 

NAURU99 

99/6/15-99/7/18 

12°S, 130°E-8°N, 167°E 

PACSF 99 

99/11/2-99/12/1 

8°S-12°N, 95-1 2 1°W 

TTWE 

91/11/21-91/12/13 

0°N, 140°W 
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Table 2. Comparison of daily LHF, surface wind speed (U), surface air specific humidity (Q a ), and 
SST of GSSTF2 with those of all nine field experiments and five tropical experiments with large 
samples. The mean is ship-observed values averaged over collocated days, positive bias indicates 
larger GSSTF2, SD error is std dev of differences, and r is correlation coefficient. Units are W m' 2 
for LHF, m s' 1 for U, g kg' 1 for Q a and °C for SST. 


Source 

Days 

Variable 

Mean 

Bias 

SD error 

Daily Monthly 

r 

all 

167 

LHF 

99.8 

0.8 

35.8 

6.5 

0.83 

Tropics 

134 

LHF 

93.0 

-2.6 

29.7 

5.4 

0.80 

all 

240 

U 

5.9 

0.36 

1.38 

0.25 

0.92 

Tropics 

139 

U 

4.6 

0.31 

1.07 

0.20 

0.87 

all 

240 

Qa 

15.7 

0.67 

1.23 

0.22 

0.97 

Tropics 

139 

Qa 

17.7 

1.01 

1.11 

0.20 

0.85 

all 

279 

SST 

25.7 

0.04 

0.51 

0.09 

1.00 

Tropics 

157 

SST 

28.4 

0.02 

0.30 

0.05 

0.99 
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Table 3. Regional-mean LHF, U 10m , Q s -Q IOm , Q 10m , and Q s over global oceans during 1992-93 for 
GSSTF2, as well as the differences of HO APS, NCEP, and da Silva from GSSTF2. 


Variable 

Source 

60°S-60°N 

20°N-60°N 

20°S-20°N 

20°S-60°S 


LHF 

GSSTF2 

108.2 

104.1 

122.1 

94.3 


(W m' 2 ) 

HOAPS 

-19.7 

-14.3 

-37.4 

-2.0 



NCEP 

-3.4 

+0.5 

+3.0 

-13.2 



Da Silva 

-8.5 

-5.3 

+1.3 

-21.8 


u 10m 

GSSTF2 

7.39 

7.33 

6.63 

8.25 


(ms' 1 ) 

HOAPS 

-0.26 

+0.22 

-1.13 

+0.45 



NCEP 

-0.59 

-0.27 

-1.04 

-0.26 



Da Silva 

+0.30 

+0.84 

+0.27 

+0.04 


Qs-QlOm 

GSSTF2 

4.14 

3.80 

5.10 

3.20 


(g kg- 1 ) 

HOAPS 

-0.48 

-0.43 

-0.70 

-0.23 



NCEP 

-0.51 

-0.67 

-0.09 

-0.92 



Da Silva 

-0.53 

-0.70 

-0.27 

-0.74 


QlOm 

GSSTF2 

12.07 

10.19 

16.63 

7.81 


(g kg- 1 ) 

HOAPS 

+0.83 

+0.71 

+1.07 

+0.59 



NCEP 

+0.85 

+0.95 

+0.53 

+1.14 



Da Silva 

+0.95 

+1.14 

+0.67 

+1.15 


Qs 

GSSTF2 

16.22 

14.01 

21.73 

11.01 


(g kg 1 ) 

HOAPS 

+0.35 

+0.28 

+0.37 

+0.36 



NCEP 

+0.33 

+0.29 

+0.44 

+0.22 



Da Silva 

+0.41 

+0.44 

+0.44 

+0.41 
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Table 4. Regional-mean std dev of differences/temporal cross correlation of monthly LHF, U 10m , 
Q s — Qjom and Q 10m between GSSTF2 and each of HO APS, NCEP, and da Silva over global oceans 
during 1992-93. 


Variable 

Source 

60°S-60°N 

20°N-60°N 

20°S-20°N 

20°S-60°S 

LHF 

HOAPS 

16.8/0.81 

18.9/0.87 

16.9/0.83 

15.6/0.77 

(W m' 2 ) 

NCEP 

18.3/0.78 

17.4/0.88 

23.0/0.69 

13.4/0.84 


Da Silva 

25.7/0.51 

25.1/0.73 

27.0/0.54 

24.4/0.35 

u 10m 

HOAPS 

0.67/0.85 

0.79/0.87 

0.61/0.87 

0.67/0.81 

(ms 1 ) 

NCEP 

0.63/0.84 

0.59/0.91 

0.71/0.80 

0.56/0.86 


Da Silva 

1.01/0.55 

0.98/0.71 

0.93/0.64 

1.12/0.35 

Qs~QlOm 

HOAPS 

0.52/0.77 

0.57/0.81 

0.55/0.81 

0.45/0.70 

(g kg' 1 ) 

NCEP 

0.71/0.50 

0.69/0.70 

0.93/0.27 

0.47/0.66 


Da Silva 

0.75/0.34 

0.78/0.59 

0.80/0.31 

0.67/0.24 

QlOm 

HOAPS 

0.38/0.97 

0.48/0.98 

0.37/0.97 

0.35/0.96 

(g kg' 1 ) 

NCEP 

0.70/0.85 

0.67/0.96 

0.92/0.71 

0.46/0.95 


Da Silva 

0.83/0.81 

0.82/0.94 

0.88/0.78 

0.77/0.79 



29 


TABLE CAPTIONS 

Table 1. Times and locations of nine field experiments conducted by the NOAA/ETL research ships. 

Table 2. Comparison of daily LHF, surface wind speed (U), surface air specific humidity (Q a ), and 
SST of GSSTF2 with those of all nine field experiments and five tropica] experiments with 
large samples. The mean is ship-observed values averaged over collocated days, positive 
bias indicates larger GSSTF2, SD error is std dev of differences, and r is correlation 
coefficient. Units are W rri 2 for LHF, m s' 1 for U, g kg’ 1 for Q a , and °C for SST. 

Table 3. Regional-mean LHF, U 10m , Q s -Q 10m , Q 10m , and Q s over global oceans during 1992-93 for 
GSSTF2, as well as the differences of HO APS, NCEP, and da Silva from GSSTF2. 

Table 4. Regional-mean std dev of differences/temporal cross correlation of monthly LHF, U 10m , 
Q s — Qiom and Q 10m between GSSTF2 and each of HOAPS, NCEP, and da Silva over global 
oceans during 1992-93. 
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FIGURE CAPTIONS 

Fig. 1. GSSTF2 daily (a) latent heat fluxes, (b) surface winds, and (c) surface air specific 
humidity versus those of nine field experiments. C is for CO ARE, F for FASTEX, and X 
for other experiments. 

Fig. 2. The 10-m wind speed averaged over 1992-93 for (a) GSSTF2, and differences of (b) 
HO APS, (c) NCEP, and (d) da Silva from GSSTF2. Only collocated monthly valid data 
from all four datasets are used. 

Fig.3. Standard deviations of differences (left) and temporal cross correlation (right) of monthly 

10-m wind speeds between GSSTF2 and each of (a, b) HOAPS, (c, d) NCEP, and (e, f) 
da Silva during 1992-93. 

Fig. 4. Same as Fig. 2, except for 10-m specific humidity. 

Fig. 5. Same as Fig. 2, except for sea-lOm humidity differences. 

Fig. 6. Same as Fig. 3, except for sea-lOm humidity differences. 

Fig. 7. Same as Fig. 2, except for latent heat fluxes. 

Fig. 8. Same as Fig. 3, except for latent heat fluxes. 
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